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Abstract 

This paper reconsiders finite variable reductions of the universal 
Whitham hierarchy of genus zero in the perspective of dispersionless 
Hirota equations. In the case of one-variable reduction, dispersionless 
Hirota equations turn out to be a powerful tool for understanding the 
mechanism of reduction. All relevant equations describing the reduc- 
tion (Lowner-type equations and diagonal hydrodynamic equations) 
can be thereby derived and justified in a unified manner. The case of 
multi-variable reductions is not so straightforward. Nevertheless, the 
reduction procedure can be formulated in a general form, and justi- 
fied with the aid of dispersionless Hirota equations. As an application, 
previous results of Guil, Manas and Martinez Alonso are reconfirmed 
in this formulation. 
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1 Introduction 



The Lowner equation pQ is a differential equation that describes a family of 
deformations (with parameter A) of univalent conformal maps g\ : D — ► B\ 
from a fixed disk D C CP 1 to a simply connected domain B\ = B\T\ C CP 1 
with a slit formed by a continuously growing arc I\ on a fixed curve T. g\ 
is the inverse of a univalent conformal map f\\B\—>D whose existence is 
ensured by Riemann's mapping theorem. For technical reasons, we assume 
that D is centered at oo, B contains oo in its interior, and g\{oo) = oo. To 
write the Lowner equation, let us express the maps g\ and f\ in terms of 
coordinates as z = g(p,\) and p = f(z,X). The Lowner equation thereby 
reads 



for f(z, A) and 



df(z, A) «(A) + /(z,A)0#A) 

<9A Jl2: ' j K (A) -/(^,A) <9A 



gg(p, ^) = P + «(A) gg(p, A) g0(A) 
<9A p — k(A) 9p 9A 



for g(p, A). (As one can easily see, these two equations are equivalent.) k(A) 
and 0(A) are auxiliary functions that are determined by (and, conversely, 
determine) the deformations of the maps. k(A) is called a "driving force", 
and 0(A) is related to the behavior of g(p, A) and f(z, A) at infinity: 

g( P) \) = e m P + 0(l) (p-oo), 
/( Z ,A) = e- m z + 0(l) (z^oo). 

A few variant of the Lowner equation are also known. Firstly, one can 
consider a family of univalent conformal maps g\ to a domain B\ with 
multiple slits of arcs on several fixed curves Ti, . . . ,T M . Since the arcs 
can grow independently, this family depends on M-dimensional parameters 
A = (Ai, . . . , Am)- The univalent homomorphic functions g(p, A) and f(z, A) 
representing g\ and its inverse fx satisfy a system of Lowner-like equations 

df(z, A) _ Kj (\) + f(z,X) d<j>(\) 

d\j nZ ' - f(z, A) dXj ' 

dg{p,\) p + Kj(\) dg(p, X) d(f>(\) . 

aAj p — Kj(A) aAj 

with respect to A/s. Secondly, one can choose the upper half plane H rather 
than the disk, and consider univalent conformal maps g\ from H to a simply 
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connected domain B\ with oo on its boundary. The holomorphic functions 
g(p, A) and f(z, A) representing these maps satisfy a pair of differential equa- 
tions called the "chordal Lowner equation" (see below). This equation, too, 
have multi-slit analogues. 

Remarkably, these Lowner-type equations also emerge in finite variable 
(or finite field) reductions of various dispersionless integrable systems. This 
fact was first discovered by Gibbons and Tsarev in the case of the Benney 
equations [21 E]. g(p, A) shows up therein as a generating function g(p) of 
the Benney moments, and the equations of motion of those moments are 
converted to an evolution equation for g(p). This evolution equation can be 
identified with the Lax equation of the second flow of the dispersionless KP 
hierarchy. Gibbons and Tsarev observed that a general M-variable reduction 
of this system can be obtained in the following form: 

1. g(p) depends on the spacetime variable (x,t) of the Benney equations 
via M reduced variables A = (Ai, . . . , Am) as 

g(p) = g(p,x), x = \(x,t). 



2. g(p, A) satisfies the chordal Lowner equations 

dgjp, A) = 1 dgjp, A) da{\) 
d\j p — Uj(X) dp d\j 



j = l,...,M, 



of domains B\ with M slits. The inverse function f(z,X) of g(p, A) 
satisfies the dual equations 

df(z, A) 1 da(X) 

dXj Uj(X) - f(z, A) dXj ' 3 ' 

Uj(X) and a (A) are auxiliary functions that determine the reduction. 
Uj(Xys are called driving forces, and a(A) is related to the asymptotic 
behavior of g(p, A) and f(z, A) at infinity. 

3. Aj's satisfy a system of "diagonal" hydrodynamic equations 

OX; ^ s dXj . 

■^ = »(A)af, , = i,...,m. 

Thus Xj's are "Riemann invariants" of the reduced system. The "char- 
acteristic speeds" Xj{X) turn out to satisfy a set of conditions, which 
enables one to solve these equations by Tsarev's generalized hodograph 
method ^j. 
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This result has been generalized to many other cases El El [HI El QUI ED] • 
Recently this issue is also studied from the point of view of Hamiltonian 
structures [T2] . 

In this paper, we address the problem of finite variable reductions of 
the universal Whitham hierarchy of genus zero [13] ■ This issue was studied 
by Guil, Manas and Martinez Alonso [5], and the relevance of Lowner-type 
equations was already recognized therein. We reconsider this issue in a new 
perspective based on the dispersionless Hirota equations. As demonstrated in 
the case of the dispersionless KP and Toda hierarchies [H] , the dispersionless 
Hirota equations can be a powerful tool for studying finite variable reductions 
of dispersionless integrable systems. The goal of this paper is to generalize 
this observation to the case of the universal Whitham hierarchy. 

This paper is organized as follows. Section 2 is a brief review of the univer- 
sal Whitham hierarchy of genus zero. Sections 3 and 4 deal with one-variable 
reduction. In Section 3, we derive Lowner-type equations and diagonal hy- 
drodynamic equations as necessary conditions for one-variable reduction. In 
Section 4, we show that the hydrodynamic equations can be solved by the 
hodograph method, and lastly confirm that the Lowner-type equations and 
the hydrodynamic equations are sufficient to determine a solution of the uni- 
versal Whitham hierarchy. Section 5, 6 and 7 are devoted to multi-variable 
reductions. In Section 5, we examine "rational reduction" (or "algebraic or- 
bits") of the universal Whitham hierarchy as a prototype of multi- variable 
reductions. Following the lines illustrated therein, we formulate a general 
form of multi- variable reduction in Section 6, and reconsider the results of 
Guil et al. [9J on the basis of this formulation. Section 8 is a summary of 
our results. 

2 Universal Whitham hierarchy of genus zero 

This section is a collection of basic notions and results on the universal 
Whitham hierarchy of genus zero picked out from the literature [TBI [141 [T5l 

Eg. 

2.1 Lax equations 

We consider the universal Whitham hierarchy of genus zero with N + 1 
marked points on a Riemann sphere with coordinate p [13j . One of the 
marked points are fixed to p = oo; the others p — q 1 , . . . , are part of the 
dynamical variables of the hierarchy. The other dynamical variables are the 
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coefficients of the Laurent series 



z 



3=1 

r, 



z a (p) = — h y^u aj (p - q a y 1 (a = 1, . . . , N), 

p-q* jr[ 

which are assumed to converge in domains T>q and T> a (or asymptotic expan- 
sion of holomorphic functions at the boundary points p = oo and p = q a of 
such domains). These Laurent series are dispersionless analogues of the Lax 
operators in dispersive integrable hierarchies. 

In this setup, the universal Whitham hierarchy has N + 1 series of time 
variables t an (n = 1, 2, . . . for a = and n — 0, 1, . . . for a — 1, ... , N). Let 
us use t to denote these time variables collectively, namely, 

t = (^01; ^02; • • • 5 ^10) til, ... , tNO, £/Vl> • • •)• 

The lowest variables £oi> tio, ■ ■ ■ , tm i n eac h series t an , a = 0, 1, . . . , N, play 
a special role. In particular, i i should be interpreted as a spatial variable 
rather than a time variable. In a sense, tio, . . . , two, too, are a kind of spatial 
variables, which originate in charge variables of an N + 1-component charged 
fermion system [16]. For convenience, we introduce the auxiliary variable 

N 
a=l 

(This stems from the constraint of total charge being zero.) We also use the 
abbreviation 

dan 9/ dt an . 

for the derivatives in the time variables. 

Time evolution of the system is defined by the dispersionless Lax equa- 
tions [T3] 

danZpip) = {n an (p), z p (p)} {a,/3 = 0,l,...,N) (2) 
with respect to the two-dimensional Poisson bracket 

, . , df dg df dg 
dp dt i dt i dp 
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on the (p,t i) space. Thus p is a conjugate variable of £ ij i n other words, 
classical limit of the "momentum" 9 i i n a one-dimensional quantum me- 
chanical system. Qon(p)' s an d ^an(p)' s for n > 1 are polynomial in p and 

(p - ^r 1 , 



^On(p) = P n + a n2p n 2 H h aOnn, 

n «n(p) = T——TZ + 7Z I W=l + • • • + 



and given by the singular part of Laurent expansion of zo(p) n and z a (p) r 
(including the constant term for the former): 



z (p) n = n 0n (p) + o( P l ) (p^oo 

z a (p) n = n an (p)+0(l) (p^q a ). 



(4) 



The first few of them read 

fioi(p)=P, fi 02 (» = p 2 + 2w 02 , 

noi(p) = -^, 0^) = -^^ + ^^, .... 

fi Q0 (p)' s are exceptional and given by logarithmic functions: 

Q a o(p) = - log(p - g a )- (5) 

2.2 ^-functions and Hamilton- Jacob i equations 

We now introduce a set of new variables to extend the foregoing Lax equations 
to a larger system. 

The first set of new variables are "conjugate variables" of zp(p)'s [THIS], 
which play the role of the Orlov-Schulman operators in the present setup. 
As a consequence of the Lax equations, f2 an (p)'s satisfy the dispersionless 
Zakharov-Shabat equations [T3] 



dp n tt am {p) - d am Q[3 n (p) + {VLam{p), Qf3n(p)} = («, (3 = 0, 1, . . . , N) . (6) 

They can be packed into the equation 

uj A uj = (7) 

for the closed 2-form 



TV oo 

UJ = 

n=l a=l n=0 



^2 d ^0n(p) A dt 0n + ^2 ^2 dtt an (p) A dt an , 
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where "d" stands for exterior differentiation with respect to both t and p. 
By Darboux's theorem, u can be expressed in a canonical form with two 
"Darboux variables". One can choose z@(p) as one of the Darboux variable; 
let (p(p) denote the conjugate variable: 

u = dz p (p)AdCf,(p) (J3 = 0,1,...,N). (8) 

This equation implies that C/3G°)' S satisfy Lax equations of the same form 

danUp) = {n*n(p), Cp(p)} 0, 1, N) (9) 

as zp(pYs, along with the canonical Poisson relations 

{zp(p),(p(p)} = l (P = 0,1,...,N). (10) 

Having introduced C/3(p)' s ? we can now introduce a second set of new 
variables Sp(p), (5 = 0, 1, . . . , N, which are called "S-functions" [T31 El El US]. 
To this end, we rewrite (jSj) as 

d(6 + Cp(p)dz p (p)) = 0, 

where 

oo N oo 



n=l a=l n=0 

The S'-function Sp(p) is defined as a potential of 6 + (/3(p)dzp(p): 

e + Cp(p)dz p (p) = dSp(p). (ii) 

Though one can find these S"-functions as Laurent series of p (for (3 = 0) or 
oi {p — q a )~~ l (for f3 = 1, . . . , N), it is more convenient to consider them as 
Laurent series in zp{p). In such an expression, they can be expanded as 

5d(p) = ^2 ^nZoipT + ho log Z (p) - ^ V 0n, 

1 1 ^ 

n=l n=l /^2 N \ 

°° °° Zfl(p)~ n 

5/j(p) = ^ tfinZfiipT + tpo log Z/j(p) + 0/3 - ^ Vp n , 

n=l n=l 

where the coefficients v 0n , vp n and 4>p are functions of i. We can rewrite this 
expression as 

Soip) = S Q (z (p)), 5/j (p) = Sp(z p (p)) 
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by introducing the functions 

oo oo 

s o( z ) = ^2 tonzTl + to ° log z ~ £ — v ° n 

n=l n=l 71 (13) 

oo oo 

Sp{ z ) = ^2 t f 3nZ;n + log z + ~ — v > 3n 

n=l n=l 

of z. These Spfy's, too, are called "S'-functions". 

In a sense, the second set of S'-functions Sp(z) are more fundamental. 
They satisfy the Hamilton- Jacobi equations [HU E] 

d an Sp(z) = Q an (d 01 Sp(z)) (a = 0,l,...,iV), (14) 

which are quasiclassical limit of the (scalar-valued) auxiliary linear equations 
of the underlying multi-component KP hierarchy [16]. Moreover, the Voi- 
der ivative 

Pp(z) = d 01 Sp(z) 
gives the inverse function of z — Zp(p), namely, 

ZpMz)) = z, PpMp)) = P- (15) 

We can derive a simple relation among qp,rp and <fip,Vpi from these re- 
marks. The foregoing Laurent expansion of the S'-functions implies that po(z) 
and pp(z) (P = 1, . . . , N) behave as 

p (z) = z + 0(z~ l ), pp(z) = d 01 (f)o + doi^i-aT 1 + 0(z~ 2 ) 

as z — > oo. Since the inverse function p = pp{z) of 

z = Zf,(p) = -^- + 0{l) 
P-Q/3 

should coincide with this expression of pp(z), we find that 

qp = -<9 O 10/3, rp = -doxvp!. (16) 

On the other hand, the Hamilton- Jacobi equations for t a0 take such a form 
as 

d a0 Sp(z) = - \og(d iSp(z) - q a ) = - \og{pp(z) - q a ), 
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which can be solved for pp{z) as 

pp{z) = q a + e- a ^ z) . 
Letting a — (3 and recalling the Laurent expansion of Sp(z), we find that 

p p (z) = q p + e-^z- 1 + 0{z- 2 ). 
This gives another expression of rp\ 

r/3 = e - da0rpa . (17) 

2.3 F-function and Hirota equations 

We now introduce the F-function T = F{t) (logarithm of the quasiclassical 
r function [15] ) as a solution of the equations 

daoF = ~(p a + log(-l), a = 1, . . . , N, 

13=1 

where log(— 1) is understood to be equal to, say, ni, though the choice of 
the branch is irrelevant in the final result. (This definition of the F-function 
[To] is slightly different from that of Manas et al. [TT1 [To] , but this does not 
affect the main part of results.) This strange factor is related to the signature 
factors € a j3 that we shall encounter below. 

All relevant quantities of the hierarchy can be expressed in terms of the 
F-function. In particular, the S"-functions have a very compact expression: 



S (Z) = ton*" + too log Z - D Q {z)F, 

71=1 

oo 

S a (z) = ^2 t an z n + t a0 log z + <fr a - F> Q (z)F, 



(19) 



n=l 



where 



00 z -n 00 z -n 

D o( z ) = ^ d on, D a (z) = Y] d an . 

n=l n=l 

Consequently, the p-functions, too, can be neatly expressed as 

Po(z) = z - d 01 D (z)J r , 

Pa{z) = -d l(D a (z) + dao)?. 



(20) 
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The ^-functions do not have such a simple expression, but one can determine 
the Laurent coefficients u an order by order. For example, let us derive an 
expression of M02, q a ,r a - The first few terms of the foregoing expression of 
the p-functions read 

p (z) = z-d 2 m Fz- l + 0{z- 2 ), 

p a (z) = -d m d aG F - d m d al Fz~ l + 0(z~ 2 ). 

On the other hand, as the inverse function of 



zoip) = P + U02P 1 + 0(p 2 ), 



Za(p) = — 2L- + 0(1), 

Po(z) and p a {z) should have Laurent expansion of the form 



po(z) = z - U02Z 1 + 0(z 2 ' 



Pa{z) = q a + r a z 1 + 0(z 2 ). 

Thus we find that 

«02 = doi-? 7 ; Qa = -dmdaoF, r a = -d m d al F. (21) 

Note that this expression of q a and r a can also be derived from (fTE]) . 

When the universal Whitham hierarchy is viewed as dispersionless (or 
quasiclassical) limit of the (charged) multi-component KP hierarchy [16J , the 
F-function is identified with the limit of logarithm of the tau function. One 
can thereby derive the "dispersionless Hirota equations" (or, more appropri- 
ately, "dispersionless differential Fay identities") 

p D Q {z)b {w)T _ 1 _ doi(A)(z) -Dq(w))F 
e 1 , 

z — w 

ze b o(z )b a{w) r = z _ doi (f) o{z) _ b a {w))T, 

p b a (z)b a { w) F _ zwd 01 (D a (z) - D a {w))F 

z — w 

e ^ e b a(z )b,( w )F = _ dm (D a {z) - b p {w))F (a + (5) 
for a, /3 — 1, . . . , N, where 

D {z) = { D ^ z ) ( a = °)' e = / +1 

a{Z> \D a (z) + d a0 (a^O), €aP 1-1 (a>(3). 



(22) 
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The signature factor e a p stems from a fermionic representation of the tau 
functions of the multi-component KP hierarchy [TTl [T5] . 

These dispersionless Hirota equations are equivalent to the universal Whitham 
hierarchy itself [16J. This is a generalization of the results known for the dis- 
persionless KP and Toda hierarchies [T9l 1201 |21| [22] . The notion of Faber 
polynomials and Grunsky coefficients, which originate in complex analysis, 
play a crucial role here. 

2.4 Faber polynomials and Grunsky coefficients 

Faber polynomials and Grunsky coefficients are hidden in the foregoing setup 
of the universal Whitham hierarchy. 

First of all, Q 0n (z) and Q an (z) are nothing but the Faber polynomials 
of po(z) and p a (z), respectively [IE]. Namely, these "polynomials" (in p 
and in {p — Qa)" 1 , respectively) are characterized by the following generating 
functions: 



Po(z -q ST z nf\ 

log = - > fton(g), 

z — ' n 



Z 1 (23) 

log = - > — n m ? . 

q-q a ^ n 

Moreover, differentiating these generating functions by q yields the generating 
functions 

oo 

1 \ - Z 



n=l 

oo 

1 v Z 



(24) 



q-Pa{z) q-q a ^ n 



for the derivatives 

0n(g) = c>g ' = dq ' 

It will be more suggestive (and convenient) to rewrite the second equation 
of flMD as 

oo _ n 

r^=-£V^-^- (25) 



n 



This type of identities are called "kernel formulas" in the literature 
(because the left hand side may be thought of as a Cauchy kernel), and shown 
to be useful in some applications of dispersionless Hirota equations such as 
"associativity equations" (2TJJ, El] . 

Grunsky coefficients show up when we use (1201) to rewrite the right hand 
side of the dispersionless Hirota equations ff2"2"]) as 

d Q1 (D (z) - bo{w))F 

1 = po(z) - p {w), 

z — w 

z - d 01 (D (z) - D a {w))T = po(z) - p a (w), 
zwd 01 (D a (z) - b a {w))F zw{p a {z) - p a (w)) 



z — w z — w 

-d 01 (b a (z) - Dp{w))J : = p a (z) - Pp{w), 

and consider the logarithm of both hand sides of each equation. (122j) thus 
turn into the following equations 

b (z)b (w)F = log 
b (z)b (w)F = log 

b a (z)b a (w)F = log 







-p {w) 




z - 


- w 


Po{ 


*)- 








• 

Z 


zw 


(Pa 


(z) - Pai 






z — w 


Pai 


» 


-pp(w) 



(26) 



As one can show using the definition of the Faber polynomials ()23l) . these 
equations are, actually, a generating functional representation of the Hamilton- 
Jacobi equations (|14P ; this fact lies in the heart of the aforementioned equiv- 
alence of the dispersionless Hirota equations and the universal Whitham 
hierarchy. 

The right hand side of ( 1261) are nothing but generating functions of the 
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(generalized) Grunsky coefficients b am p n of the p-functions: 
log Po(z)-PoM = _ A z . mw _ nhQ 

7 — in < J 



OmOn , 

m,n=l 
oo oo 



Po(z)-p a (w) -,„ 

log £ — ^ — = - 2^ 2^ -2 w 



Oman j 

m=l n=0 

oo 



Zw( Pa (z) - Pa (w)) y Z ~ m W -% : 

z — w ^ 



(27) 



m,n=0 

p a (z)-pp(w) ^ _ m 

log^^ PV_^ = _ ^ z m w ( a ^ 

m,n=0 

In particular, the z~ x terms of the first and second equations yield the iden- 
tities 

oo oo 
W - p {w) = ^ w ~ nb 010n, ~Pa{w) = - W n b i an . (28) 
n=l n=0 

One can thus recover the p-functions from the Grunsky coefficients. 

( I26p show that all second derivatives of the F-function are given by the 
Grunsky coefficients as 

damdpnF = -b am (3n («, = 0, 1, . . . , N), (29) 

where we have introduce the rescaled derivatives 



^d an (n^O), 
d a0 (n = 0). 



Conversely, one can use ( 1291) as defining equations of F. In that case, the 
Grunsky coefficients have to satisfy a set of integrability conditions: 

■ynf3m ■ 

(30) 

Once those integrability conditions are shown to be satisfied, one can obtain 
the F-function as a solution of these equations. Moreover, this simulta- 
neously ensures that the dispersionless Hirota equations are also satisfied, 
because ( 12U]) are, after all, the dispersionless Hirota equations in disguise. 
This method for showing the existence of the F-function was first developed 
for the case of the dispersionless KP and Toda hierarchies [TT]. We shall 
apply this method to reductions of the universal Whitham hierarchy. 
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3 Lowner-type equations in one-variable re- 
duction 



In the setup of one-variable reduction, the dynamical variables u an (n = 
0, 1, . . ., u a0 = r a ) and q a are assumed to depend on t via a single function 
A = A(t) as 

u an = u an (X(t)), q a = q a (X(t)). 
Consequently, z a (p) and p a (z) are thought of as functions of p and X(t): 

z a (p) = z a (jp, A(t)), p a (z) = p a (z, A(i)). 
The goal of this section is to derive the following 

Theorem 1. If u an and q a are functions of a single variable X = X(t), 
then there is a function U = U(X) of X such that p a (z) = p a (z,X) and 
Zaip) = z a (p, X), a = 0,1, . . . , N , satisfy the Lowner-type equations 

dp a (z) 1 du Q2 



(31) 
(32) 



dX U-p a (z) dX ' 

dz a (p) z' a (p) du 02 
dX p-U dX' 

where z' a (p) denotes the p- derivative 

z' (p) = 
a[P) dp ■ 

Moreover, X = A(i) satisfies the hydrodynamic equations 

9 an X = Xan{X)d m X (33) 

with characteristic speeds 

Xan (X) = tf m {U). 

This theorem is a generalization of the one- variable reduction of the KP 
and Toda hierarchies [TTJ. Let us explain some implications: 



1. In a quite general situation as the theorem assumes, we can say nothing 
about the shape of the domain of p a (z) (equivalently, the range of z a (p)) 
and the position of U therein. In particular, it is a priori not evident 
whether the domain of p a (z) is a slit domain like the one assumed in 
the original setup of the Lowner equation. Since the reduction defined 
by (l3"Tj) and (!3"2l seems to be meaningful in such a general situation (cf. 
various solutions of Lowner-type solutions known in the literature p3, 
El [7J El El Ell E] ) , we do not restrict our consideration to slit domains. 
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2. All equations of (13T1) and f[3"2"j) have a common driving force U. In 
particular, (13"2"1) mean that 2 a (p)'s satisfy an identical linear differential 
equation of the form 

d_ 1 du 02 d \ ,^_ n 
p-U d\ d P J Za[P) ~ 

By the classical theory of characteristics, this implies that z a (p)'s are 
mutually functionally related. Namely, there are functions f a (z) of one 
variable z such that 

*<*(?) = /«(*>(?)), « = l,...,iV. (34) 

This agrees with the setup of reductions by Guil et al. [9] , though they 
assume the special form 

/«(*) = , (35) 

where Cry S cLFG constants. 



3. (|33p determine the time evolution of the reduced dynamical variable 
A = A(i). The whole hierarchy is thus reduced to the hydrodynamic 
equations ( |33l) with Riemann invariant A and characteristic speeds 
X«n(A). These hydrodynamic equations can be solved by the hodo- 
graph method (see Theorem EJ). 

4. It is not obvious from the subsequent proof that the converse of the 
statement of the theorem also holds. Therefore one has to show, sepa- 
rately, that (I3T1) (or (1321 ) and (|33|) lead to a solution of the universal 
Whitham hierarchy. We shall do it with the aid of the dispersionless 
Hirota equations (see Theorem [H]). 

3.1 Proof of Theorem [D step 1 

Differentiating both hand sides of the first equation of fl26l) by t i yields the 
equation 

o n /x n / <9oiPo(^) - d 0l p (w) 

dQ l DQ{z)D Q {w)J r = — — . 

Po(z) -po(w) 

Note that D (z) = D (z). Let us rewrite this equation as 

doiPo(z) - d 01 p (w) 



Po{z) -Po{w) 



d 01 D (z)D (w)F 

doiPo(z) d 01 p (w) 



D (w)p (z) D (z)p (w) 
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We have used (1201) to rewrite the denominator d iD (z)D (w)J-' as 

d iD {z)D (w) T = -D (w)po(z) = -D (z)p (w). 
By the chain rule, we can express the derivatives on the right hand side as 

dp (z) dp (w) 
ooiPo{z) = 7^ OqiA, a iPo(w) = — d iA, 

D (w)po{z) = — D (w)X, D (z)p {w) = D (z)X. 
Hence last equation reduces to 

Po{z) - Po{w) = - + 



D (w)X D (z)X 
or, equivalent ly, 

/ \ doiX , N c^oiA 



D (z)X y DoHA" 

Therefore both hand sides of this equation are independent of z and w. Let 
Uq = Uq(X) denote this quantity: 

Po{z) ~ T^lzjx = u °- (36) 

On the other hand, applying D Q (z) to both hand sides of the first formula 
of (12T1) and using (120]) yield the identity 

D (z)u 02 = dl l D (z)T = -d 01 p (z). 

Again by the chain rule, the derivatives on both hand sides can be expressed 
as 

D (z)u 02 = -^r-D (z)X, d 01 po{z) = — d iA 
Thus we find that 

d 01 X du 02 /dX 



D (z)X dp (z)/dX 

We can use ( 1371) to rewrite the foregoing equation ( 1361) as 

du 02 /dX 
Po{z) + r, . wm = C/b- 
9p (-2)/oA 

This implies that po(- 2; ) satisfies the Lowner-type equation 

dpo(z) _ 1 dup2 
dX U — Po(z) dX 



(37) 



(3? 
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3.2 Proof of Theorem Q} step 2 

We can repeat almost the same calculations for the second equation of ( j26l ). 
Firstly, differentiating both hand sides by toi, we obtain the equation 

a n r^n /■ Mr $>ii>a(*) - 3ni>a(w) 
a i J D a (2;) J D Q (w)J r = 



p Q (z) -p a (w) 
which can be rewritten as 

doip a (z) - d 01 p a (w) 



Pa{z) -p a {w) 



d iD a {z)D a {w)J : 
d ip a (z) d ip a (w) 



D a (z)p a (z) D a (w)p a (w) 
By the chain rule, this equation reduces to 

PM ~ a ; . = Pa{w) - 



D a (z)X D a (w)X 

This implies that both hand sides are actually independent of z and w. Thus 
we have the equation 

p a (z) - J^- = U a , (39) 

where U a = U a (X) is a function of A only 

Secondly, we can derive, from the first formula of (ED), the identity 

D a (z)u 02 = d^D^z)? 7 = -d Q1 p a (z). 

By the chain rule, this identity reduces to 

d iX du Q2 /dX 



D a (z) d Pa (z)/dX' 
391) and ( HUj) imply that p a (z) satisfies the Lowner-type equation 

dp a (z) _ 1 du 02 

dX U a — p a (z) dX 



(40) 



(41) 
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3.3 Proof of Theorem Q} step 3 

Let us examine implications of the third and fourth equation of ( !26l) . 
Differentiating the third equation by toi yields the equation 

/ x /x doiPo(z) d 01 p a (w) 

P0{Z) - Pa{W) = —z— — + 



D a (w)po(z) D (z)p a (w) 
By the chain rule, this equation reduces to 

du 02 /dX du 02 /dX 

Po{Z) ~ Pa{w) = - - + 



do (z)/dX dp a {w)/dX 

By ( 1381) and ( HT1) . we can rewrite the two terms on the right hand side as 

du 02 /d\ du 02 /d\ 

- U - po{z), ^ — . . , =U a - p a {w). 



dp (z)/d\ ' dp a (w)/d\ 

Thus we find that 

U = U a , a = l,...,N. (42) 

namely, (138jl and ( T4T1) actually have an identical driving force U. 

The same conclusion follows from the fourth equation of ( 1261) . Namely, 
differentiating this equation by £oi eventually leads to the identities U a = Up. 

3.4 Proof of Theorem [D step4 

We can derive the evolution equations ( 1331) as follows. 
Let us rewrite ( 1371) as 

nn , dp (z)/dX d 01 X 
Dq{z)X = — 7^t-5oiA - 



du 02 /dX U - p (z) 

and use the identity 

1 



U-po{z) ^ n 



that can be obtained from the kernel formula ( 1241) by letting q = U. The 
outcome is the equation 



D (z)X = J2—^on(U)d 01 X. 



n 

n=l 
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This is a generating functional form of the equations of f[3"3"j) for a = 0, 
n = l,2,.... 

In exactly the same way, using ( |40l) and the kernel formula for p a (z), we 
can derive the equations 

/ 00 -n \ 

D a (z)X = (D a (z) + d a0 )X = IJ2 —V<*n(U) + V a0 (U)j d 01 X. 

They are a generating functional form of the equations of f[3"3"j) for a = 
1, . . . , N and n = 0, 1, 



4 Solutions from one-variable reduction 
4.1 Hodograph method 



The hydrodynamic equations (13311 can be solved by the hodograph method. 
The hodograph method is extremely simplified in this case, because there is 
only one variable A. 

Theorem 2. Let F(X) be an arbitrary function of X, and X = X(t) a function 
that satisfies the hodograph equation 

co N 00 

^t nX0n(A) + ^^t an Xan(A) =F{\). (43) 
n=l a=l n=0 

Further assume that the regularity condition 

V% dX0n(A) ^Xan(A) ^ gF(A) ^ 

2^ 0n ax + 2^ 2^ an ax ^ ^a~ (44) 

n=l a=l ri=0 

/ioWs for X = X(t). Then X = X(t) satisfies the hydrodynamic equations ( fffgj) . 

Proof We differentiating both hand sides of the hodograph equation by 
t an . By the chain rule, this yields the equations 

f\\ , f V% ^Xom(A) <9x^ m (A)\ dF{X) . 

Xanl^J + I / v tOm r / r /3m 7^ I Can^ — — — Can^, 

\m=l 1 (3=1 m=0 / 

hence 

v - f dF(A) A ^Xom(A) ^n(A)^ 

\ m=l ' j3=l m=0 J 

19 



In particular, letting a = and n = 0, we have the equation 

_ / dF(X) ^ gxom(A) v>v>. dxp m (\) \ a . 

\ m=l /3=1 m=0 / 

If we multiply the last equation by XanW and subtract it from the previous 
one, the outcome is the equation 



OO r, / , n N OO 



9F(A) _ v - dxo m (X) x^x^. dxpmj X) 

ox ^ 0m d\ 2^2^ l v™ dx 

m=l /3=1 m=0 



/3m " 



X {d an X - Xan(X)d 01 X). 

By the regularity condition, we can drop the prefactor of d an X — Xcm(A)<9 iA 
and obtain the hydrodynamic equations ( 1331) . □ 



4.2 Existence of F-function 

In the last section, we derived the Lowner-type equations (EH), ( 1321) and 
the hydrodynamic equations ([3"B"]) . but we have not confirmed the converse, 
namely, whether these equations ensure that u an = u a (X(t)) and q a = 
Qa{X(t)) give a solution of the universal Whitham hierarchy. 

We now prove that the converse is also true, following the idea presented 
in the end of Section 2.4. This is also a generalization of the result in the 
case of the dispersionless KP and To da hierarchies [11] . 

Theorem 3. The integrability conditions fW\) of are satisfied in the 
foregoing setup of one-variable reduction. The F -function T = ^(i) thus 
defined by 12$) gives a solution of the dispersionless Hirota equations |Hj]. 



Proof Let us first illustrate the calculations in the case where a = (3 = 7 = 
0. We substitute z = z\ and w = Z2 in the first generating function of ( 1271) . 
and apply D (z 3 ) = Dq{z 3 ) to both hand sides. This yields the generating 
function 

D (z 3 )(p ( Zl ) - p (z 2 )) ^ , _ m A 

7 \ / s — — / , z l z 2 Z 3 °0n°0l0m 

Po(Zl) ~ Po(Z2) 

of donboiom's. On the other hand, recall that the hydrodynamic equations 
(1331) have a generating functional representation (cf. step 4 of the proof of 
Theorem [1]). In the case of a = 0, it reads 

D (z)X = —— -. 

Po{z) - U 
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Using this equation, the Lowner-type equation for po(z) and the chain rule, 
we can rewrite the left hand side of the foregoing generating function of 
donboiom's as 

P0O1) ~Po(z 2 ) 

D (z 3 )X ( dpo(zi) dp (z 2 ) 



Po(zi) -Po(z 2 ) V 9X dX 

d 01 X (I 1 \ <9w 02 



(Po(zz) - U)(pq(zi) - pq(z 2 )) \U - pq(zi) U-p (z 2 )J dX 

doiX du 02 

(U - p ( Zl ))(U - Po (z 2 ))(U - Po (z 3 )) dX ' 

Since this quantity is symmetric in z\ and z 3 , we have the functional identity 

A)(z3)(Po(*i) ~Po(z 2 )) _ Do(zi){po(z 3 ) -Po(z 2 )) 
PoW) -Po{z 2 ) Po(z 3 ) -Pa(z 2 ) 

This implies the identities 

donboiom = d ib 

of the coefficients, which are exactly the integrability conditions ( 1301) for 

a = (3 = 7 = 0. 

In much the same way, staring with one of the generating functions of 
(127)) and applying D 7 (z) to it, we obtain the generating functions 



^ Z\ Z2 z^ d^ n b 



Lm,n 



of derivatives of the general Grunsky coefficients. Since the hydrodynamic 
equations (1331) have the generating functional form 

D a (z)X 



p a (z) - U' 

we can rewrite the foregoing generating functions as 

D-r(z3){Pa(zi) -Pp{z 2 )) doiX du 02 

Pa(zi) - Pp(z 2 ) (U -p a {z x )){U -p p {z 2 )){U -p 7 (zg)) dX ' 

This implies the functional identity 

fl 7 (*3)(Pa(zi) -Pfjjzi)) = D a (z 1 )(p 1 (z 3 ) -pp(z 2 )) 
Pa(zi) ~ Pp(z 2 ) p 7 (z 3 ) -Pl3{z 2 ) 
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hence the integrability conditions flHUj) as expected. More precisely, we have 
to be careful about the difference of the four types of generating functions in 
( 1271) . but this does not affect the final conclusion. □ 
Let us mention that this proof reveals an interesting feature of the inte- 
grability conditions ( 1301) . Namely, since the derivatives of the Grunsky coeffi- 
cients are nothing but the third derivatives of the F-function, the functional 
identities in the proof imply that these third derivatives have a generating 
function of the form 

l,m,n 

dpi A 

(U-pMXU-ppteMU-p^za)) d\ ■ 

This result is suggestive from the point of view of associativity equations 
[201 El], because the third derivatives are fundamental quantities therein. 
We shall see that a similar result holds in multi- variable reductions as well. 

5 Rational reductions (algebraic orbits) 

Unfortunately, it seems difficult to extend the foregoing method for the one- 
variable reduction to multi-variable reductions. In the case of multi- variable 
reductions, we shall start from Lowner-type equations rather than derive 
them. The problem is how to find a correct form of Lowner-type equations. 
Though an answer to this question is presented in the work of Guil et al. [9] , 
we dare to take a different (heuristic) route that leads to the same answer. 

Our strategy is to examine rational reductions of the universal Whitham 
hierarchy (or "algebraic orbits" in the terminology of Krichever [TS]) as a 
prototype of general multi-variable reductions. This class of reductions cover, 
for example, the Zakharov reduction of the Benney equations [25] , reductions 
of the dispersionless KP hierarchy related to 2D topological field theories 
[261 [271 128] . a hydrodynamic reduction of the Boyer-Finley equation [29], 
etc. 

The work of Ferapontov, Korotkin and Shramchenko [2H] is particularly 
suggestive, because their method is exactly based on Lowner-type equations 
and some related equations. They indeed used those equations to apply 
Tsarev's generalized hodograph method [1]. On the other hand, since they 
do not use Lax equations, one cannot readily see how to generalize their 
results to higher flows of an underlying hierarchy. Therefore it is essential to 
understand their method in the perspective of Lax equations. 



du 



02 



(45) 
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Bearing these issues in mind, let us briefly look into rational reductions of 
the universal Whitham hierarchy. Let us mention that the following consid- 
eration is more or less parallel to the approach that Gibbons et al. followed 
in the case of the Benney equations [21 [3J, [5] . 



Setup of rational reduction In a rational reduction, we assume that 
there is a rational function 

kg N k a 



E( P ) = /° + E a °nP ko ~ n + E E 



n=2 



i i (P - Qa) n 
a=l n=l 



with poles at p = oo, q±, . . . , gjv such that zo(p) and z a (p), a = 1, . . . , N, are 
given by Laurent expansion of fractional powers of E(p) as 

z o (p) = (Laurent expansion at p = oo) 

z a (p) = E(p) 1 ^ ka (Laurent expansion at p = q a ). 

Actually, E(p) can have extra (dynamical) poles other than oo and g a 's. This 
is indeed the case for, e.g., the Zakharov reductions of the Benney equations 
and the hydrodynamic reduction of the Boyer-Finley equations. In those 
reductions, E(p) takes such a form as 

A/-1 



£(p)=p+E^V' 



whereas the Benney equations and the Boyer-Finley equations are embedded 
into the one- and two-point universal Whitham hierarchies (in other words, 
the dispersionless KP and Toda hierarchy). 
In this setup, the critical points 

p = p v E'{ Pj )=0, j = l,...,M, 

and the critical values 



X j = E(p j ), j = l,...,M, 

of E(p) play the role of driving forces and Riemann invariants. More precisely, 
if E(p) is sufficiently general, the critical values A/s can be used as full 
parameters (or "moduli") of E(p). Thus E{p) is understood to be a function 
E(p, A) of p and A = (Ai, . . . , Xm)- 
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Hydrodynamic equations for Xj The Lax equations ([2]) for the ^-functions 
reduce to the Lax equations 



d an E{p) = {Q an ,E(p)} (a — 0,1, ... , N) 
for E(p). Written more explicitly, these equations read 

d an E( P ) = n' an (p)d iE(p) - E'( P )d 01 n an (p). 

Letting p = pj in this equation yields the equation 



(46) 



d an E{p) 



p=Pj 



n' an ( Pj )d 01 E( P ) 



p=pj 



On the other hand, by the chain rule, differentiating Xj = E(j>j) by t an gives 
the identity 



p=Pj 



d an Xj = E'(pj)d an pj + d an E(p) = d an E(p) 

p=Pj 

Thus we obtain the diagonal hydrodynamic equations 

d an Xj = XanWdoxXj (a = 0, 1, ... , N) (47) 
with characteristic speed 

Note that p/s are now understood to be algebraic functions pj = Pj(X) 
defined by the equation E'{p) = 0. Thus (T4T1) may be thought of as a closed 
evolutionary system for A. As we shall show in a more general case, these 
reduced equations can be solved by the generalized hodograph method. 

Lowner-type equations for E(p) We now consider E(p) to be a function 
of p and A = A(t), and use the chain rule to rewrite the Lax equations ( 1461 ). 
Both hand sides the Lax equations can be thereby expressed as 



t ttq \^ 9E(p) A dEjp) 

LHS = 2^ -q^- dan A,- = 2^ ~dx- Q ^{P. 

j=l 3 



and 



^)d iXj 



24 



where we have used ( H7I) as well. Thus the Lax equations reduce to 

E ((^(P) - <UPi))^ - ^(P)^^) 9 oiA, = 0. 
Consequently, if = £?(p, A) satisfies the equation 



gives a solution of the 

A=A(t) 



then for any solution A = \(t) of ( )47j) . -E(p) 
Lax equation. 

Let us examine (118]) in more detail. When a = and n = 1, this equation 
is a trivial identity. The lowest nontrivial one is the case where a = and 
n = 2. In this case, since fi 02 (j») = + 2w 2, reduces to 

/ ^9E(p) du 02 n 

or, equivalently, 

ggg) = E^_du^ (49) 
9Aj p — pj d\j ' 

This is exactly a Lowner-type equation. We can readily derive the equations 

dzgjp) = z f a {p) du 02 , 5Q . 
d\j p — pj d\j 

for the Zq,(j>)'s, because the z-functions are given by fractional powers of 
E{p). The p-functions, in turn, satisfy the dual equations 

dp a (z) _ 1 du 02 
^Aj Pj — p a {z) d\j 

As regards the case for n > 2, fHHl) turn out to be redundant, namely, 
automatically satisfied if ( H9l) is satisfied. This fact can be explained in a 
more general form; we shall return to this issue in Section 6.2. 



6 Multi-variable reductions 

In view of the foregoing interpretation of rational reductions, it is now rather 
straightforward to find a correct formulation of general multi-variable reduc- 
tions. In this section, we present this formulation and its implications. After 
all, this reduction procedure is nothing but the "diagonal reduction" in the 
sense of Guil et al. [9j. We, however, attempt to reformulate it along the 
lines that we have pursued in the case of one-variable reduction. 
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6.1 Lowner-type equations for p a (z) and z a (p) 

In an M-variable reduction, the fundamental dynamical variables u an (n = 
0, 1, . . ., u a o = r a ) and q a are assumed to be functions of M-dimensional 
reduced dynamical variables A = (Ai, . . . , Am)- The reduced dynamical vari- 
ables, in turn, depend on t as A = X(t) and eventually satisfy a set of 
diagonal hydrodynamic evolution equations. 

Accordingly, thep-functionspo^z) and the ^-functions z a (p), a — 0, 1, . . . , N, 
are functions of p and A, 

z a (p) = z a (p, A), p a (z) = p a (z, A). 

We assume that these functions satisfy the Lowner-type equations 

dp a (z) 1 du 02 . 



d\j Uj — p a (z) 8\j 

or the dual equations 

dz a (p) __ z' a (p) du 02 
d\j p — Uj d\j ' 



M, (52) 



j = l,...,M, (53) 



for a given set of driving forces Uj = Uj (A) and the auxiliary function m 2 — 
Uq2 (A). Unlike the one- variable case, these auxiliary functions have to satisfy 
a set of integrability conditions. We shall present these conditions later on 
when we consider the hodograph method. 

As in the case of one-variable reduction, the dual equations ( 1531) imply 
that the z-functions are functionally related with each other by functions 
f a (z) of one variable z as flH^|) shows. Though Guil et al. [5] further assumed 
the special form ( 13*51 . the following consideration is not limited to that case. 

As a consequence of these Lowner-type equations, z a (p)'s turn out to 
have critical points at p — Uj, namely, 

z' a (U ) = 0. (54) 

This is an immediate consequence of the structure of ([531 . Since the left 
hand side has no singularity at p = Uj, the pole of l/(p — Uj) at p — Uj has 
to be canceled by a zero of z' a (p). This, however, does not imply that the 
critical values z a (Uj) coincide with A/s. This is in accord with the general 
fact that the choice of Riemann invariants is not unique but allows large 
arbitrariness. A standard way will be to choose the critical values zo(Uj) 
of Zq(z) as Aj's; the critical values zJJJj) of the other z-functions are then 
functionally related to z (Uj) as (IM|) implies. 
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6.2 Hydrodynamic equations for Aj 

In view of the results on one-variable and rational reductions, it seems plau- 
sible that the reduced variables A = A(t) satisfy diagonal hydrodynamic 
equations of the form 

WdoiXj (55) 

with characteristic speeds 

The relation between the Lax equations (EJ) and these equations, however, 
is more delicate than in the case of rational reductions, because A/s in the 
present setup are not assumed to be given by the critical values of the z- 
functions. Therefore we cannot derive (|55p by simply letting p = Uj in the 
Lax equations. Nevertheless, we can confirm that (|55|) are correct equations 
to be satisfied by A = A(t). 

Theorem 4. If / Tjgj) (or §E2\)) an d are satisfied, then the Lax equations 
ID) are also satisfied. 

Although this theorem can also be deduced from some other results that 
we shall show later on (Theorem [6] and Theorem [7]) , we dare to present a 
direct proof here. The outline of the proof is parallel to the case of rational 
reductions. A clue is the following 

Lemma 1. If p a (z) ? s satisfy §EE) , then fl an {p) ! s satisfy the identities 
dVL an (p) VL' an (p) - Kn(Uj) du 02 



d\j p — Uj d\j ' 



(56) 



Proof Let us first consider the case where a = 0. We differentiate the 
generating function (I2"3j) of f2 0n (p)'s by Aj and use (15*21) . This yields the 
identity 

z- n dVL 0n (p) 1 dp (z) 



-E 



^ n dXj Po(z) — V 9Xj 

1 du 02 

(p - Po (z))(Uj - po(z)) dXj ' 

On the other hand, the kernel formula ( 1241) implies another identity of the 
form 



DC 



~—n i 



n - ■ p-Po(z) Uj -po(z) 



Uj-p 



(p-po(z))(Uj -Pq{z))' 
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By comparing these identities, we find the identity 



M z- n dn 0n ( P ) = ^ *-» n> 0n (p) - n' a0 (Uj 

n=l J n=l 



^ n dXj ^ n p — Uj 



which implies that ( 1561) holds for a = 0. In much the same way using the 
second kernel formula (|23|) . we can derive (l5l*"j) for a = 1, . . . , N. □ 



Proof of Theorem [4] The reasoning in the case of rational reductions also 
works in this case with slightest modification. Firstly, using (I55p and the 
chain rule, we can rewrite the Lax equations ([2]) as 

E ((«Uri - - 4&> )^^) *^ = »■ 

On the other hand, combining the equations (f5"3l for ^(p) with the identity 
of the lemma yields the equations 

dz p (p) _ z'pip) dn an ( P ) 



dx j K^-KniUj) ox 
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which implies that the last equations are indeed satisfied. □ 
Let us note that the last equations (|57j) amount to (pfBl) in the previous 
section. As mentioned therein, fHHj) contains the Lowner-type equations (j49]) 
as a special case with a = and n = 2. In the present setup, ff5*3"]) is a special 
case of floTl) with a = and n = 2. By the way, Lemma - ]] (or its proof) says 
that ( 1571) is a consequence of ( 1531) . This explains why ( )48j) are "redundant" 
in the setup of the last section. 



6.3 Hodograph method 

The hydrodynamic equations (l55l) can be solved by Tsarev's hodograph 
method [3]. As it turns out below, this is a straightforward generalization of 
the framework developed by Gibbons and Tsarev [21 E] for reductions of the 
Benney equations. 

In the mult i- variable case, (15*2"]) and (15*31 have to satisfy a set of integra- 
bility conditions. As regards ( l52l) . the integrability conditions can be derived 
by eliminating the A-derivatives of p a {z) from the identity 

d ( 1 du 02 \ d ( 1 du 02 \ 
dXj \U k -p a (z) dX k ) dX k \Uj -p a (z) dXj J ' 
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After some algebra, these conditions reduce to the equations 

dU k 1 du 02 



dX 3 



Uj — U k d\j 



d 2 u 02 _ 2 du 02 du m 



(5f 



which take exactly the same form as the equations derived by Gibbons and 
Tsarev in the case of the Benney equations. The same equations can be 
derived from the integrability conditions 







du 02 d d 



1 du 02 d 



dXj p — Uj dXj dp dX k p — U k dX k dp 







(59) 



of the dual equations (joHj) as well. 

An important consequence of these equations is the following 

Lemma 2. The characteristic speeds Xanj = XanjW of K55\) satisfy the 
equations 



dx 



ank 



dX 3 



(.Xanj Xank)Vj 



where 



du 2 



(IJ. _ Uk )2 QX 3 ■ 



(60) 



(61) 



Proof Let us first consider the case where a = 0. Substituting q — U k in 
the kernel formula (|24p . yields the identity 



-E 

71=1 



z- n _ 1 

n U k -p (z) 
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We now differentiate both hand sides by Xj. Using (J52]) and ([58]) . we can 
rewrite the outcome as 

z~ n dxonk 1 f dU k dp (z 



^ // 0\, (V,-i>,(-))' 2 \ox, ox, 



I f I 1 \ 9m 2 



(C4 - j9 (^)) 2 \Uj - C/ fe Uj-poiz); 8X 3 

1 du 02 

(U k - p (z))(Uj - Po (z))(U, - U k ) dXj 
1 1 \ 1 du 2 



U, - p (z) U k - p (z) J {U 5 - U k f dX 3 



^2 n ^anj 



Xank / 



n=l 



1 du 02 

(Uj - u k f ax j • 



This shows that (l60l) are indeed satisfied for a = 0. We can confirm (l60l) for 
a — 1 , . . . , N in the same way, now using the second kernel formula (125j) . □ 
We can now formulate the hodograph method for ( 1551) as follows. 



Theorem 5. Lei = i^-(A) ; j = 1, . . . , M, 6e a set of functions of X that 
satisfy the equations 

^ = (Fj - F k )V jk , (62) 
and X = X(t) a solution of the hodograph equations 

oo N oo 

5^*0nX0nj(A) +^5^t<mXcmj(A) = Fj(X), j = 1, . . . , N. (63) 
n=l a=l ra=0 

Further assume that the regularity conditions 

2^ ton ^x~ + ^ ^ tan dx 3 ^^a~' J " 1 '---' iV ' (64) 

n=l J a=l n=0 J J 

hold for X = X(t). Then X = X(t) satisfies the hydrodynamic equations (E3J). 

Proof We differentiate both hand sides of (I6"3"j) by t an . By the chain rule, 
this yields the equations 

m / oo dx ■ N dx \ M dF- 

fc=l \m=l fc /3=1 m=0 fc / fc=l fe 
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hence 

fe=l \ " m=l " 0=1 m=0 fc / 

Let us examine the quantity inside the parenthesis on the right hand side. If 
k 7^ j, we can use (160]) and ( 1621) to rewrite this quantity as 

r, „ OO r, N OO „ 



K m=l re /3=1 m=0 



/3m " 

AT 00 



{Fk — Fj)Vkj — ^0m(X0mfc — X0mj)Vkj — t(3ra{X(3mk — X/3mj)Vkj 

m=l /9=1 m=0 

00 AT 00 \ 

i^fe - ^ tomXomk tfimXfimk - (k replaced with j) J V kj 

m=l (3=1 m=0 / 



which vanishes by flboj) . Thus the last equation simplifies as 

Xanj [dXj ^ 0m d\, ^^ f3m d\ j anr 

\ J m=l J (3=1 m=Q ■> / 

Moreover, if a = and n = 2, this equation reduces to 

' r>p OO N OO \ 



J m=l J /3=1 m=0 J 



If we multiply the last equation by Xanj an d subtract it from the previous 
one, the outcome is the equation 

' r\ 7-1 OO o N OO o 



E^-EE 



tflrn ^ X 



m=l /3=1 m=0 ■> 

x (d an \j — XanjdaOl^j)- 

By the regularity condition, we can drop the prefactor of Xan^j — Xcmjd a m\j 
and obtain the hydrodynamic equations ( 1551) . □ 
We are thus eventually left with the problem of finding F/s that satisfy 
). Such functions are given by contour integrals of the form 

A r dp G a (z a {p)) 
j = 2^t ™ ( n -n,v > (65) 



a=0 



Ca 2m (p-Utf' 
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where G a (z) is an arbitrary holomorphic function of one variable z defined 
on the range of the map p t— > z a (p), and C a is a closed curve (or cycle) in 
the domain of z a (p). It is not difficult to show that these F/s satisfy ( !62l) as 
a consequence of ( |53i) . Note that the characteristic speeds Xanj themselves 
have such a contour integral representation as 



where the path of integral is a small circle encircling the point p = oo (a = 0) 
or p = q a (a = 1, . . . , N). 

6.4 Existence of F-function 

Theorem [3] and its proof can be generalized to the present setup without 
substantial modifications. 

Theorem 6. The integrability conditions (3U\) of are satisfied in the 
foregoing setup of multi-variable reduction. The F -function T = ^(i) thus 
defined by / T^Pj) gives a solution of the dispersionless Hirota equations ((H)). 

Proof We can proceed just as in the proof of Theorem [3J By applying 
D*y(z) to the generating functions fl2T|) . we obtain the generating functions 



of derivatives of the Grunsky coefficients. By the kernel formulas f l2^|) and 
( |25l) . the hydrodynamic equations ( 1551) can be cast into the generating func- 
tional form 



We can thereby rewrite the foregoing generating function as 




(66) 




£>7(*3)(Pa(*i) -Pp{z 2 )) 

Pa(zi) ~Pp(z 2 ) 




M 
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This implies the functional identity 

Pa(zi) ~ Pp(z 2 ) Pj(z 3 ) ~ Pl3{z 2 ) 



The integrability conditions flHOl) follow from this identity immediately. The 
rest of the statement of the theorem is a consequence of the comments in the 
end of Section 2.4. □ 
As a byproduct of this proof, we obtain the following generalization of 
151) to multi-variable reductions: 



^ Z\ Z2 d a idp m dry n J~ 



l,m,n 

m n , „ (67) 



-E 



doiXj du 



02 



i (Uj - p a (zi))(Uj - p/s^iUj - p 7 (z 3 )) 



7 S'-functions in multi-variable reductions 

As an application of the foregoing formulation of mult i- variable reductions, 
we now reconsider the construction of S-functions by Guil et al. [§]. As it 
turns out below, some part of their construction can be made more trans- 
parent with the aid of the identities ( |56|) of Lemmadl 

The construction of S"-functions by Guil. et al. is based on hodograph 
solutions of the hydrodynamic equations ( 1551) . Given such a solution along 
with the p-functions satisfying (!52*|) . they construct the S'-function in such a 
form as 

00 N oo 

S P(P) = *0n^0n(p) + ^2 ^n^an(p) + S/3- (p) (68) 
n=l a=l n=0 

or, equivalent ly, 

oo N oo 

^C 2 ) = 5j*On^On(P/3(^)) + t an Q an (pp(z)) + Sp-{jpp{z)) (69) 

n=l o=l n=0 

where <Sg_ (p) are required to satisfy the equations 

/ d 1 du Q2 d 

\ d\j p — Uj d\j dp 



) VM = (70) 
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Since the differential operators on the left hand side are commutative, see 
(JSH]), the integrability conditions of these inhomogeneous Lowner equations 
are given by 







dX 



j P 







1 du 02 d 
- Uj d\j dp 
1 du 02 d 



du 



02 



p-U k dX k 



du 



02 



dX k p-U k dX k dp J \p- Uj dXj 

By straightforward calculations, one can see that these conditions are equiv- 
alent to the equations for F/s. Thus the existence of a solution to (ITU]) 
is ensured in the setup of the hodograph solution. With these definitions, 
the main part of the results of Guil et al. can be stated as follows: 

Theorem 7 (Guil, Manas, Martinez Alonso |9J). The S -functions Sp(z), 
(3 = 0, 1, . . . , N, satisfy the Hamilton- Jacobi equations [14\ ). 

Let us prove this theorem using our tools. Firstly, by the chain rule, the 
tan-derivative of Sp(z) can be expanded as 



M 



j=l 7,m 

+£(Vm*)>^ 



dXj 



+ 



+ 



dXj 



P=Pg(z) 
r) \ 



dan Xj 



(71) 



where we have used the abbreviated notation 

CO 



N co 
^ ^ 1 



■ym 



m=l 



urn + 

7=1 m=0 



Using the Lowner-like equations and the identity (IoT)l) of Lemmedl we 
can calculate the quantity in the first parenthesis on the right hand side of 



tt'(Pf3(z)) 



dpp{z ) t dVL im (p) 



+ 



P=Pp{z) 



%mW^) du 02 % m {pp{z)) - fyJJJj) du 02 

U-pp(z) dXj pp(z)-Uj dXj 

%m( U j) dUp 2 

pp{z)- Uj dXj ' 
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As regards the quantity in the second parenthesis, we use the equations (170]) 
satisfied by (z) as 



Since each term of the sum on the right hand side vanishes by the hodograph 
equations (1031 . we have the equation 

d an S p {z) = Vt an (pp(z)). 

In the case where a = and n = 1, this equation reduces to 



by which we can eliminate pp{z) from the last equation and obtain the 
Hamilton- Jacobi equations ( THj) . 

This completes the proof of the theorem. Note that using the identities 
(I56p makes the proof shorter and more understandable than the original proof 
of Guil et al. 

8 Conclusion 

We have thus seen that Lowner-type equations play a fundamental role in 
finite variable reductions of the universal Whitham hierarchy of genus zero. 
The status of dispersionless Hirota equations therein is more subtle. As 
regards the one-variable reduction, the dispersionless Hirota equations (J22J) 
are certainly a clue. The generating functional form (1261) of these equations 
enabled us to derive the Lowner-type equations (13~TI) and (|32j) directly from 
the assumption that all dynamical variables are functions of a single reduced 




Pp{z) - Uj d\j ' 



Consequently, (ITT]) turns into an equation of the form 




doiSp(z) = P/3{z), 
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variable A. Unfortunately, this method does not work for multi-variable 
reductions. Therefore we were forced to start from (rather than derive) the 
Lowner-type equations (l52l) and ( !53l) and to confirm that they are a correct 
set of reduction conditions (Theorem HI). For both one- variable and the 
multi-variable reductions, however, we could eventually justify the reduction 
procedure in a unified way, namely, by proving that the defining equations 
( l2"9"j) of the F-function is integrable (Theorems El and E]). This is a place 
where the dispersionless Hirota equations play a truly fundamental role. 

Viewed from a technical point of view, another clue of our method is the 
use of the generating functions (1231) and ( 1271) and various identities derived 
therefrom. Not only being closely related to the dispersionless Hirota equa- 
tions themselves, these generating functions turned out to be also extremely 
useful in many aspects of finite variable reductions. 

Our next target will be, naturally, the cases of nonzero genera jT3J . It 
will be rather easy to derive dispersionless Hirota equations for those cases, 
as partly argued by Krichever et al. in a different setup |30j. A main prob- 
lem is to find a correct form of Lowner-type equations. We expect to find a 
prototype in differential geometry of Hurwitz spaces [3TJ E2] and associated 
Whitham-type hierarchies [T3l 127], because they are generalizations of the ra- 
tional reductions that we considered as a prototype of general multi- variable 
reductions for the genus zero case. Presumably, we should start with the 
genus one case, for which an explicit description of Hurwitz spaces are avail- 
able in the literature [3T1 ESI EH EH] along with a candidate of Lowner-type 
equations [36] . 
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